EPJ manuscript No. 

(will be inserted by the editor) 



Axial charges of excited nucleons from Cl-fermions 

T. Maurer\ T. Burch\ L. Ya. Glozman^, C. B. Lang^, D. Mohler^ and A. Schafcri ^ 

^ Institut fiir Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, Germany 
^ Institut fiir Physik, FB Theoretische Physik, Universitat Graz, A-8010 Graz, Austria 
3 TRIUMF, 4004 Wesbrook MaU Vancouver, BC V6T 2A3, Canada 

the date of receipt and acceptance should be inserted later 

Abstract. We report lattice QCD results on the axial charges of ground and excited nucleon states of both 
parities. This is the first study of these quantities with approximately chiral (CI) fermions. Two energy 
levels in the range of the negative parity resonances TV* (1535) and A'^* (1650) are observed and we determine 
the axial charge for both. We obtain a small axial charge for one of them, which is consistent with the 
chiral symmetry restoration in this state as well as with the small axial charge of the A'^*(1535) predicted 
within the quark model. This result agrees with the findings of Takahashi et al. [1] obtained with Wilson 
quarks which violate chiral symmetry for finite lattice spacing. At the same time for the other observed 
negative parity state we obtain a large axial charge, that is close to the axial charge of the nucleon. This 
is in disagreement both with the quark model prediction as well as with the chiral restoration but allows 
for an interpretation as an s-wave tt N state. 
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1 Introduction 

With the advent of high luminosity accelerators hadron 
spectroscopy became again part of the forefront of physics. 
New states have been observed, e.g., the X,Y,Z bosons by 
Babar, Belle and Cleo [HIS] which suggest that states exist 
which have a more complicated lowest-Fock-state parton 
structure than qq (and correspondingly qqq). One exten- 
sively discussed possibility is the existence of tetraquarks. 
One of the issues is where to draw the line between tetra- 
quarks and meson molecules [4l|5], or pentaquarks and 
baryon-meson molecules for that matter. In this context 
the two negative parity nucleon excitations iV*(1535) and 
A^*(1650) are of special interest as it is often assumed that 
they are mixtures of a genuine three quark state and a pen- 
taquark state, which can also be thought of as a molecule 
[51[71l51[51ITUj . So, a better understanding their properties 
is a topical task. Another one is to improve on hadron 
spectroscopy in general, and thus to be better able to de- 
cide which hadrons do not fit into standard qq and qqq 
phenomenology. One of the major order schemes for the 
latter is chiral symmetry. 

If chirality were a good symmetry of hadron physics, all 
hadrons would occur in parity pairs, which is obviously not 
the case. However, one can hope that remnants of chiral 
symmetry remain. Therefore, it is a much discussed ques- 
tion whether hadron resonances show some statistically 
relevant degree of parity pairing. The answer seems to be 
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affirmative, see, e.g., [TT], but not indisputable. Therefore, 
the identifications and measurement of quantities which 
are sensitive to this issue are important. 

In refs. p^lT^ it was proposed that the observed par- 
ity doubling for higher excitations can be explained via an 
effective restoration of chiral symmetry and the formation 
of approximate parity-chiral multiplets. Furthermore, the 
axial charge of each parity doublet state has to be zero 
for an exact symmetry. The hypothesis of an effective chi- 
ral symmetry restoration can therefore be tested by cal- 
culating the axial charges of members of parity doublet 
resonances [T3]. 

There exists, in fact, circumstantial evidence for the 
hypothesis of effective chiral restoration. The pion decay 
rate of nucleon resonances with approximately restored 
chiral symmetry should be suppressed and, in fact, the 
decay of probable members of parity doublets was found 
to be suppressed by a factor of ten pij. 

For this and many similar problems in hadron physics 
it is hoped that lattice QCD will provide crucial clari- 
fications. In particular, the variational analysis method 
15,16J has helped considerably to study excited states 
on the lattice 1171118,19,20.21,22,23,24 . However, finding 
the correct interpretation of the results of such simulations 
is always a challenge. 

Ideally, given large enough propagation distance in Eu- 
clidean time, high enough statistics and a complete set 
of interpolators, the observed energy levels in the prop- 
agator analysis should be well defined. In reality there 
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always remains a certain dependence on the chosen set 
of interpolators. Typically, only those states which have 
a large overlap with one or more of the set become sta- 
tistically significant. Using just three-quark interpolators 
for the baryons, it is unclear whether one can observe a 
noticeable overlap with baryon-meson states, even in sim- 
ulations with dynamical quarks. As an example for such 
weak coupling one may refer to the p-meson, where only 
explicit inclusion of p-wave TTTr-interpolators lead to cor- 
responding energy levels in the analysis The stabil- 
ity of the eigenvectors over several time slices helps to 
identify the eigenstate; however, in the absence of explicit 
baryon-meson interpolators, one cannot hope to discern 
the physical content of the state. 

Indeed, the expected lowest energy for the s-wave TriV 
channel overlaps with the region of the two lowest nega- 
tive parity 1/2" states, observed on the lattice by practi- 
cally all groups. Consequently, one does not know a priori 
whether the observed two lowest states in dynamical lat- 
tice calculations represent the resonances A^*(1535) and 
A^*(1650) or, in fact, one of the resonances and the tt ~ N 
system in the s-wave. One needs additional dynamical in- 
formation about these states to resolve the issue. Axial 
charges of the observed states could help to answer this 
question. 

A study similar to our's was already done some time 
ago by Takahashi et al. [55] who indeed found that one 
of the first negative parity nucleon resonances has a very 
small axial charge. However, these authors used Wilson 
fermions which violate chiral symmetry for any finite a, 
which obviously is a cause of concern for this kind of in- 
vestigation. To clarify the situation we present here a sim- 
ilar study, performed, however, with approximately chiral 
fermions, so-called CI (chirally improved) fermions. The 
motivation for the development of this fcrmion action \27\ 
[2H] is that implementing chirality exactly on the lattice 
introduces unavoidable non-localities which make simula- 
tions very expensive. (Such simulations became possible 
recently, see [25], although only on small lattices.) 

The outline of the paper is as follows: In section [2] we 
will discuss gA and its calculation on the lattice. In section 
[3]the variational method used will be presented and lattice 
details are given in section [4] We present and interpret the 
numerical results in section [s] and finally summarize our 
conclusions. 



2 The axial charge qa 

In this section we summarize well known facts about the 
axial charge which provides a background for what will be 
discussed further on. 

The axial charge gA of the ground state nucleon, or 
more precisely the ratio gA/gv is a central property of 
this baryon. It describes its coupling strength to the weak 
interaction and has been deduced, e.g., from neutron /3 
decay [5U] : 

g^/gy = 1.2670 ± 0.0030 



The neutron /3 decay involves four form factors (see |31j): 
vector gy, tensor g^, axial gA and pseudo-scalar gp: 

{p\V+\n) = Up [j^gviq'^) ~ qxc^xt^griq'^)] Un (1) 
(pI^m I"-) = [7/^755^(9^) - iQtilsgpiq'^)] Un (2) 

where V^^ — u'^^d, — u'^fj^j^d and is the momentum 
transfer between proton and neutron. In the limit of zero 
momentum transfer the axial and vector form 

factors dominate. Their values in this limit are called the 
vector and axial charges of the nucleon: gy = (7\/(0) and 
gA = 5a(0). 

We neglect the mass differences of up and down quarks 
(rritj — iTid ^ 0) and hence the neutron and proton mass 
difference. In this case the global chiral SU{2) x SU{2) 
symmetry (when m = m„ = rud = 0) is broken except for 
the vector SU (2) subgroup. The associated vector charge 
(7y = 1 is still conserved. The axial symmetry is explicitly 
broken by m ^ 0. The lattice calculation of axial charges 
relies on the conservation of vector symmetry. One can 
write: 

{p\V+\n) = 2{p\V^^\p) = (p|T/;b) - {p\V^\p) (3) 
{p\Al\n) = 2{p\Al\p) = {p\Al\p) - {p\Al\p), (4) 

where Vjl^' — qjfj,q, A^^^ = 97^ 759. This links the weak in- 
teraction properties to properties determined by the quark 
content of the state. The electric charge of the neutron is 
zero, which implies: 

= (nlf^n) = l{n\V-\n) - l{n\V^\n) (5) 

^ 2(n|y"|n) (nlV^ln) (6) 

Isospin symmetry implies the equivalent equation: 

2{p\V'\p) = {p\V-\p). (7) 

Since the electric charge of the proton is one, we can con- 
clude 

{p\fnp)-l{p\Vlp)-l{p\V'\p) (8) 

^{p\V-\p)~{p\V''\p)^{p\V+\n). (9) 

This implies = 1 for the nucleon. 

We now look at continuum matrix elements of states k 
of momentum p and spin s, where we normalize = — 1 
and p^ = —rn^., Generally, one defines Lorentz decompo- 
sitions 

(A:|y^('')|fc) = {k\q^^q\k) = 2v[''''''>p^ (10) 

(fc|4«)|A;) = (fc|g7M75g|fc) = ai''''h,,m, (11) 

The operator product expansion relates such matrix ele- 
ments to moments of structure functions: 

= dx [q{x) + {-irq{x)] = (12) 
Jo 

al«) = / dxx" [Aq{x) + (-l)"Z\g-(x)] = (x")^,. (13) 
Jo 
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Here, is the 0th moment of the unpolarized quark 

distribution function q{x) in a state k. It counts the num- 
ber of quarks q. For nucleons we have 

gv = v[-^-v['^ = l = {l).,^,. (14) 

Next, Oq'"''''' is the 0th moment of the polarized quark dis- 
tribution function Aq{x) and determines the average spin 
fraction carried by all quarks q. 

2gA = 4"'' ~ "■'o '' = ^Au-Ad- (15) 

If we work with states of zero momentum, we get direct 
access to these interesting quantities: 

(plK)^""'^|p) = 25ypo (16) 
{p,s\A^^''^\p,s)^2gAmpS, (17) 
In lattice normalization (p, s\p, s) = 1 one has 

{p\V^''''^\p)^gv (18) 
{p,s\A'^-'^\p,s)=gAS. (19) 

where — —1. 



3 Variational Method 

On the lattice one automatically obtains results for nucle- 
ons of opposite parity. The periodicity of the lattice allows 
propagation backwards in time, which due to T ~ T^C, 
corresponds to the propagation of antinucleons of opposite 
parity. So, for < t < Lt the two-point correlator is domi- 
nated by different parity states for small and large t values. 
For both regions one gets a superposition of ground and 
excited states. To extract specific excited states from this 
superposition is in general difficult because in Euclidean 
space-time they are exponentially suppressed with respect 
to the ground state by a factor exp{—AEt), where AE is 
the energy difference. The standard tool to deal with this 
problem is the variational method [T3], which has proven 
to be more reliable than multi-exponential fits. In this 
method one chooses a set of interpolators Oi with the 
quantum numbers of the state of interest and constructs 
a correlation matrix 

C,,{t) = (O.(t)O,(0)). (20) 

This matrix has a decomposition 

C^,{t) =Y,{mW){n\d,\Q)e-'^- . (21) 

n 

It can be shown p^, that the eigenvalues of the general- 
ized eigenvalue problem 

C{t)Mt) = ^k{t,to)C{to)Mt) (22) 

behave as 

Afe(i,to) = Cfee-(*-*°)^'= (l-t-0(e-(*-*«)^^^)), (23) 



while Afc(to,to) — 1- At fixed tg, AEk is given by 

AEk = inm{E„, -En\m^n}. (24) 

For the special case of t < 2to and a basis of N correlators 
[Sg AEk is given by 

AEk = En+i — En- (25) 

Therefore, at large time separations, each eigenvalue is 
dominated by a single state, since the difference E^^i—En 
is large. This allows for stable exponential fits to the eigen- 
values. Consequently, for large enough Euclidean time t, 
the largest eigenvalue decays with the mass of the ground 
state, the second largest eigenvalue with the mass of the 
first excited state, and so on. 

The effectiveness of this method depends crucially on 
the choice of operators. These must be sufficiently diverse 
to span a functional space large enough to have good over- 
lap with all states of interest. At the same time the basis 
cannot be chosen too large as this results in large numer- 
ical fiuctuations. 

We used a basis of six nucleon interpolators, combining 
two different spinor structures with three different levels 
of smearing (narrow, medium and wide). The two spinor 
structures are 

N^{t,p) = ^e-'P^e°''=riUa(a;,i) {u\x , t f E^d^ {x , t) 

X 

-d''{x,tfE2u''{x,t)), (26) 
where table[2]lists the chosen Dirac matrices. Interpolator 



type 


A 




Xi 


1 


C75 


X2 


75 


C 



Table 2. The structure of the two variants of nucleon inter- 
polators used. C denotes charge conjugation. 



Xi contains a scalar diquark (this scalar "diquark" in the 
local interpolator has nothing to do with the possible clus- 
tering in the physical state) , while interpolator X2 contains 
a pseudo-scalar diquark. Interpolator xi couples to the 
nucleon ground state, while the negative parity states are 
seen with interpolators xi and X2 f34llT8] . 

The three levels of smearing are indicated by a su- 
perscript /, /. More precisely, the two-point-function of a 
nucleon with momentum p, parity ±, created at time 
with quark smearings / — (Ii, 12,13) and annihilated at 
time t with quark smearings I — {Ii, I2, 13) is denoted by: 

Cii{t,p) ^ {Ni{t,p)Ni{0,p)) (27) 

In order to extract the expectation value of some op- 
erator O in a hadron state we determine the 3-point cor- 
relation matrix 

CfJr^t) ^ {B.{t)O{T)B,i0)) (28) 
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lattice 


/3 


a mo 


a/fm 




aruAwi 


mAwi/MeV 




A50 


4.70 


-0.050 


0.15032(112) 


0.7619(57) 


0.03027(8) 


39.73(10) 


200 


C64 


4.58 


-0.064 


0.15831(127) 


0.8023(64) 


0.02995(20) 


37.33(25) 


200 


C72 


4.58 


-0.072 


0.15051(115) 


0.7627(58) 


0.01728(16) 


22.65(21) 


200 


C77 


4.58 


-0.077 


0.14487(95) 


0.7342(48) 


0.01054(19) 


14.35(26) 


300 



Table 1. Parameters of our dynamical CI configurations. 



and compute ratios (see the simplest variant from |35| ) 



Rk{T,t) 



Ck{k\0\k)e-'^'' 



-tEh 



{k\0\k) 



(29) 
(30) 



to cancel the exponential factors. These ratios should be r 
independent, leading to a plateau- type behavior between 
time zero and t. The quality of this plateau is the stan- 
dard test for the precision one has obtained. Altogether 
we calculate 

± _ _ ± _ T:rP^r^y{N{t)j^y{T)Nm 

9a,v - "^AygA.y — ^ 



TTP±{N{t)N{0)) 



(31) 



with Fy = 7^ and Fj^ = 75, and the parity projectors 

l + 7t 



P+ = 



(32) 



(We choose the time direction t = 1.) 

For the different sources the corresponding matrix of 
three-point functions reads 

C^a^"'^''(T,t,p) = {Niit,p)j[{T)NiiO,p)). (33) 
with the flavor neutral current operators at time r 

^Y.4i^^^)^i'Us'i^^r). (34) 

X 

Only connected contributions are relevant to our flavor 
non-singlet matrix elements. 

4 Lattices 

We present results for four ensembles of 16'^ x 32 lattices 
and used the tadpole-improved Liischer-Weisz gauge ac- 
tion, the CI Dirac operator and np = 2 dynamical quarks, 
see table [1] (For more details on the setup cf. [T71IM] .) 

Lattice spacings have been calculated in |33[ using the 
method proposed by Sommer [571155] where the Sommer 
scale was chosen to be tq = 0.48 fm. AWI-masses extrap- 
olated to the physical point have been determined in |33j . 

We employed one step of stout smearing [35] and three 
steps of hypercubic-blocking [40 on the gauge configura- 
tions, then calculated the two- and three-point functions. 

In order to improve the overlap of our operator with 
real spatially extended states, we choose Jacobi smearing 



lattice 


identifier 


K 


N 


VRMs/a 


C77 


1 


0.223 


15 


2.30(1) 


C77 


2 


0.184 


70 


4.67(7) 


C77 


3 


0.15 


10 


0.665(2) 


C72 


1 


0.28 


7 


1.604(3) 


C72 


2 


0.1925 


37 


3.493(33) 


C72 


3 


0.4 


2 


0.8052(6) 


C64 


1 


0.28 


7 


1.603(4) 


C64 


2 


0.1918 


37 


3.463(28) 


C64 


3 


0.4 


2 


0.8052(8) 


A50 


1 


0.223 


15 


2.30(1) 


A50 


2 


0.184 


70 


4.66(8) 



Table 3. Variety of Gaussian smearings employed. Errors de- 
note standard deviations. 



[^I42| . Table [s] shows the values used, along with the 
resulting RMS radii. The capital letter always represents 
the /? value, the two numbers symbolize the value of the 
bare quark mass parameter arriQ. 

In order to reduce contamination from back-propagating 
states, we used a Dirichlet boundary condition at i^, i.e., 
we modified the gauge configuration by setting all gauge 
links from Id io to + 1 to zero. 



UtitD,x,y,z) = 



(35) 



before inverting any sources. (This was also done in [T].) 
When we calculated the lattice pion correlator we noticed 
deviations from a perfect exponential decay behavior at 
times \t — to] < 6. For this reason we chose to = 21 such 
that there is no influence from the boundary condition in 
the interval t G (0..15). We used an implementation of the 
EigCG inverter ^ to speed up to calculation of prop- 
agators for different sources on the same configuration. 
For each configuration we calculated standard propaga- 
tors and sequential propagators with all source smearings 
listed in table [3] We calculated sequential propagators for 
each T £ {2,3,4} and F £ {74,7375} to allow a deter- 
mination of matrix elements of and and hence gv 
and g^. Results for gy are only available for configuration 
ensembles C64 and C72. 



5 Numerical Results 

5.1 Nucleon Masses 

We look at positive and negative parity nucleon masses. 
From an inspection of diagonal elements of the correla- 
tion matrix Cij we find; For the positive parity nucleon, 
only interpolator xi heavily couples to the ground state 
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for all available smearing types. Interpolator X2 couples to 
excited states predominantly. For the negative parity nu- 
cleon, interpolators with xi and X2 couple to the ground 
state for all available smearings. For the negative parity 
states we observe weaker signals than for the positive par- 
ity states. We calculate the full correlation matrix dj and 
employ the variational method. 

The use of a reference matrix Cq = C(io) in the varia- 
tional method, discussed in [T^ITC] . degrades the eigen- 
value signals at the benefit of improved signal separa- 
tion for excited states. We use to = 1 for the positive 
parity states since the signal is stable enough for such a 
tradeoff. For the negative parity states we use the varia- 
tional method as well, but employ its trivial variant where 
Co = t. 

We diagonalized C{t) for every t and plotted the result- 
ing eigenvectors and eigenvalues, as well as corresponding 
effective masses. For each of the configurations we have 
full matrices and results available. 

We employed the same smearing parameters as the 
authors of ^6]. Additionally we added a third smearing 
whose RMS radius is half the size of smallest one. We 
found no effects on the masses, but the plateaus of three- 
point over two-point functions become more stable when 
we include this very narrow smearing. 

In the following we discuss only the results for the 
lowest two states of each parity because we did not obtain 
stable mass plateaus for the higher excited states. We also 
show only results for the C77 and C64 ensembles. Those 
for the other two look very similar. 

We plot results for all these subsets in figures [T] [2] and 
[3j|4j The color-coding in all of these figures is as follows: In 
the (a) and (b) figures the two lowest energy eigenstates 
are always drawn in black and red. In the (c) and (d) 
figures the coefficients for source xi and the three differ- 
ent smearings used are presented in green (narrow), black 
(medium), red (wide). Those for the source X2 are dark 
purple (narrow), blue (medium), light purple (wide). 

Unfortunately our limited computer resources allowed 
us to generate only 200 configurations for each ensemble, 
which is clearly not enough for a precise study of excited 
states. In this case, however, the precise values of gA for 
the negative parity states is not so relevant, but rather 
their approximate size. The question we want to answer is 
whether there is a negative parity state with an unusually 
small gA- Inspecting the plots Figs. [T] to |4] one should not 
only look at the mass plateaus but also at the eigenvectors 
obtained by solving the generalized eigenvalue problem. If 
these vectors are stable one has a reasonable approxima- 
tion of the eigenstate in question, even if the mass plateaus 
look marginal. With this in mind we conclude that the 
only state which is precisely extracted is the positive par- 
ity ground state. Its eigenvector composition is dominated 
by spinor structure xi with a mixture of smearing radii 1 
and 2. 

For negative parity we clearly obtain two distinct states, 
characterized by quite different eigenvectors, the masses of 
which are, however, almost degenerate with rather shaky 
mass plateaus. This is similar to our previous results, as 



ensemble 


parity 


state 


fit interval 


am 


C77 


-(- 


1 


2-11 


691f47"l 


C77 


_i_ 


2 


2-6 


1 564-(4-S') 


C77 




I 


1-5 


1.07(19) 


C77 




2 


1-5 


1.08(14) 


C72 




1 


2-11 


881f21~) 


C72 


_l_ 
i 


2 


2-6 


1 779(90"; 


C72 




1 


1-5 


1.451(79) 


C72 




2 


1-6 


1.350(63) 


C64 


+ 


1 


2-11 


1.033(6) 


C64 


+ 


2 


2-6 


1.806(18) 


C64 




1 


1-5 


1.65(2) 


C64 




2 


1-5 


1.59(2) 


A50 


+ 


1 


2-11 


0.974(18) 


A50 


+ 


2 


2-6 


1.585(35) 


A50 




1 


1-7 


1.19(16) 


A50 




2 


1-7 


1.340(85) 



Table 4. Nucleon masses measured. Fit intervals include data 
points at their boundaries. 



well as to results of other groups. Typically these two ob- 
served states are interpreted as the resonances iV*(1535) 
and A^*(1650). However, at our quark masses the S-wave 
TriV state lies in the same region. Consequently, we do not 
know whether these two observed lattice states represent 
the resonances iV*(1535) and A''* (1650), or one of these 
resonances and the 5'-wave ttN state. 

The first-excited positive parity state (the Roper) is 
not well reproduced. Its mass is much too high and its 
eigenvector composition is rather unstable. This is a prob- 
lem shared by many other lattice studies [53]. Very re- 
cent results suggest that these problems are caused by 
finite volume artefacts, which, if so, should be especially 
large for our comparably small lattices. This could indi- 
cate some special properties of the Roper wave function 
which are not yet fully understood (e.g., a strong coupling 
to more-extended Fock states). 

In order to extract the masses of the states we per- 
formed a double exponential fit to 

\{t) = Aie-'"!* + yl2e-"^* (36) 

for each set of eigenvalues, such as to correct for the higher 
state admixtures clearly visible at small t. Table|4]contains 
always the lower mass value for each such fit. 

To compare to physical masses one still has to per- 
form a chiral extrapolation. As our statistics are insuffi- 
cient to do this with reasonable precision, we simply plot 
the masses we obtained as a function of the axial Ward 
identity (AWI) quark mass in figure [5] 

5.2 Baryon Charges 

In Figs. |6] - [11] we show our results for the vector and 
axial charge for nucleons of both positive and negative 
parity. The three curves in the plots differ by the value of 
T chosen in [31] All curves should reach the same plateau 
for times t > t, which they do. The values r = 3, 4 and 5 
correspond to the black, red and green lines, respectively. 
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Fig. 1. Results for the two lowest positive parity states for ensemble C72. 
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Fig. 2. Results for the two lowest positive parity states and ensemble C64. 
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(d) Eigenvector 2, C72 



Fig. 3. Results for the two lowest negative parity states for ensemble C72. 

3 1 \ \ \ \ r- 




1 2 3 4 5 6 7 8 9 10 11 12 13 
t/a 

(a) Eigenvalues, C64 



2 3 4 5 6 7 8 9 10 11 12 13 
t/a 

(b) Effective masses, C64 




12 3 4 

(c) Eigenvector 1, C64 



5 6 7 
t/a 



9 10 11 12 13 




1 



3 4 5 6 7 8 9 10 11 12 13 
t/a 

(d) Eigenvector 2, C64 



Fig. 4. Results for the two lowest negative parity states and ensemble C64. 
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Fig. 5. Dependence of the nucleon masses on the AWI quark 
mass (ensembles C77, C72 and C64). Shown are the lowest two 
states of positive and negative parity. 

In each plot we show ratios of three-point over two- 
point functions for the state composition determined by 
the eigenvector obtained from the variational method. The 
operators chosen are 74 for the vector charge and 7375 for 
the axial vector one. We will comment on the renormal- 
ization factors below. 

In all cases presented the behavior is the same for times 
above the chosen t: the vector and axial vector charges 
form plateaus at values which are compatible for all en- 
sembles and all choices of t and thus should be physical. 
However, even when one gets a reasonable plateau, data 
are quite noisy and the list of results given in table [5] 
should be interpreted with care. Still, it seems clear that 
one of the negative parity states has indeed a very small 
axial charge. This confirms the lattice results of [T]. 

We determined the renormalization factors 

9T = Za9a, (37) 
= Zv9v. (38) 

in a separate analysis [|44j. 

As a check of the accuracy of our calculations we give 
in Table [5] also the renormalized vector charges for two of 
our ensembles. The deviation from 1, which is of the order 
0.05 should be taken as estimate for the systematic un- 
certainties of our calculation, which is substantially larger 
than the purely statistical error given in brackets. We thus 
assume that also our results for gj^ have an systematic 
uncertainty of about 0.05 in addition to the statistical un- 
certainties quoted. 

As stated above our data are insufficient for simulta- 
neous continuum and chiral extrapolations. Since we have 
data for three quark masses for ensemble C, we can at- 
tempt a linear extrapolation, see Fig.[T2|and table[6j This 
gives at least a rough feeling for the size of the effects of 
a chiral extrapolation. Table [6] gives the main results of 
our investigation. The value for for the ground state 
A^(939) comes out a bit too low, which again is a common 
problem for many lattice calculations of this quantity. It 
is known that gA is very sensitive to finite size correc- 
tions [45 | l31| . l46 1 l47 ll48ll45ll49|. l50| . which reduce its value. 



excited nucleons from Cl-fermions 
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Table 5. Renormalized axial and vector charges. Renormal- 
ization factors were taken from |44) . 



This fact is again related to larger higher-Fock-state ad- 
mixtures, which in this case may also be described in terms 
of the pion cloud [5T]. As our lattices are rather small 
{L ^ 2.5 fm) these finite size effects are a likely origin of 
that discrepancy. 

We observe two approximately degenerate negative- 
parity states in the region of the iV*(1535) and TV* (1650) 
resonances. The most important point is that we indeed 
obtain one negative parity state with a very small value of 
gA- Such a small value is consistent with approximate chi- 
ral restoration, that was initially suggested for the highly 
excited states [T^[T5l[Tl] . If the negative parity state that 
we observe is indeed a signal of the approximately restored 
chiral symmetry, then there must be its chiral partner of 
positive parity with small axial charge in the same energy 
region. This could be a possible scenario for the Roper 
resonance that has so far escaped clear identification on 
the lattice. 

The quark model predicts for the axial charges of the 
resonances TV* (1535) and TV* (1650) the values -1/9 and 
5/9, respectively [52ll53] . Consequently our small axial 
charge is also consistent with the quark model axial charge 
of TV* (1535). 

The axial charge of the other observed negative parity 
state is large, close to the nucleon's axial charge. Thus our 
result is inconsistent with result of [IJ. It is also inconsis- 
tent with both the chiral restoration scenario in this state, 
as well as with the quark model prediction. A similarity 
of the axial charge of this state with the nucleon's axial 
charge suggests the following interpretation of this state. 
The axial charge of the ttTV system in the 5-wave state 
should be close to the nucleon's axial charge since the ax- 
ial charge of the pion is 0. This hints that this state is in 
fact not a resonance, but rather a ttTV system in a relative 
S'-wave motion. Indeed, the energy of such a continuum 
state at present quark masses is the same as the ener- 
gies of both our observed negative-parity states [35] . For 
a firm conclusion on this issue one needs additional dy- 
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Fig. 6. Axial charges for ensemble C77, for explanations see text 
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Fig. 7. Axial charge for ensemble C72 



T. Maurer et al.: Axial charges of excited nucleons from Cl-fermions 
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Fig. 10. Vector charge for ensemble C72 
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Fig. 11. Vector charge for ensemble C64 
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namical observables for both negative-parity states: e.g., 
their magnetic moments. 
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parity 


state 
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Table 6. Axial charges of excited nucleons on ensemble C77. 
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6 Conclusions &. Perspectives 

We presented the first calculation for qa of negative-parity 
nucleon states using dynamical, approximately chiral (CI) 
fermions. We found a negative-parity state with small ax- 
ial charge, as predicted by the chiral symmetry restoration 
hypothesis [T^] . Such a small axial charge is also consistent 
with the axial charge of the iV*(1535) predicted by the 
quark model. This observation confirms the lattice results 
obtained by Takahashi et al. pj with Wilson fermions, 
which suggests that the chiral symmetry violation inher- 
ent for Wilson fermions has no significant effect for that 
calculation. From our calculation one cannot identify this 
state with either the iV*(1535) or the iV*(1650). Our re- 
sult for the ground state nucleon gA is rather close to the 
physical value. We attribute the remaining discrepancy to 
finite volume effects. 

The second observed negative-parity state has practi- 
cally the same axial charge as the nucleon. This suggests 
that this state may not be a resonance, but a ttN system 
in the S-wave of relative motion. Further studies of this 
interesting issue are required for a firm conclusion. 

As this was a pioneering study with dynamical CI 
fermions, it can be improved in many ways by future work. 
Most obviously, our analysis was based on an very small 
number of configurations and very few ensembles which 
strongly limited our possibilities to control the chiral, con- 
tinuum and infinite-volume extrapolations. Substantially 
more statistics would clearly help. However, of equal im- 
portance is the optimal choice of sources for such studies. 
The results of our analysis suggest that there is room for 
improvement. Actually, this is not just a technical issue. 
From the different overlap with various sources one can 
try to deduce at least qualitative information about the 
structure of the wavefunctions of hadron states, informa- 
tion which one can hardly get in any other way. More 
specifically, our results suggest that the wavefunction of 
the negative-parity state with very small axial charge is 
compact because it has large overlap with a very narrow 
source. This is consistent with our interpretation that the 
other negative-parity state represents a tt N state. To ob- 
tain more insight, it would also be interesting to repeat 
this analysis with other sources, in particular with sources 
which contain derivatives. Derivative sources have been 
shown to improve lattice determinations of excited-meson 
properties; see [541155]. 
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